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In this paper we investigate the perturbations by a massless Dirac spinor of a 
Born-Infeld black hole. The decay rates of the spinor field which is given by the quasi- 
normal mode frequencies are computed using the sixth order WKB approximation. 
The behavior of the quasi normal modes with the non-linear parameter, temperature 
and charge of the black hole are analyzed in detail. We also compare the decay rates 
i/-) ' of the Born-Infeld black hole with its linear counter-part Reissner-Nordstrom black 

hole. 
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When black holes are perturbed, the resulting oscillations lead to quasinormal modes 
(QNM). Such modes are defined as the solutions to the related wave equations char- 
acterized by purely ingoing waves at the horizon. In addition, one has to impose 
boundary conditions on the solutions at the asymptotic regions as well. In asymptot- 
ically flat space-times, the second boundary condition is for the solution to be purely 
outgoing at spatial infinity. Once these boundary conditions are imposed, the re- 
sulting frequencies become complex and discrete. The frequencies of QNM's depend 
on the parameters of the black holes such as mass, charge and angular momentum 
and are independent of initial perturbations. If the radiation due to QNM modes 
are detected in the future by gravitational wave detectors, it would be a clear way of 
identifying the possible charges of black holes. There are extensive studies of QNM's 
in various black-hole backgrounds in the literature. See the review by Kokkotas et. 
al. PQ for more information. 
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An important aspect of QNM studies have been related to the conjecture relating 
anti-de Sitter(AdS) and conformal field theory (CFT) [2J. It is conjectured that the 
imaginary part of the QNM's, which gives the time scale to decay the black hole 
perturbations, corresponds to the time scale of the CFT on the boundary to reach 
thermal equilibrium. There are many work on AdS black holes in four and higher 
dimensions on this subject and we will refer a few among them in [3] [I] [5] [6] . 

Another reason for intensive research on computing the QNM's has been due to 
Loop quantum gravity. The claims relating the real part of the high overtones of 
QNM's to entropy in loop quantum gravity has led to many works such as [7j [8] [9] 
[IO][n][12][13][Tl][T5][^[nj[l8]. 

QNM's of gravitational and scalar perturbations are extensively studied. In con- 
trast, the number of works on spinor perturbations around black holes are smaller. 
QNM's of spinors of the Schwarzschild black hole (with and without the cosmological 
constant) were studied by numerous authors. Some of them are the works by Zhi- 
denko [IS], Cho [2D], Shu et.al. [21], Jing [22] [23]. The Reissner-Nordstrom black 
hole ( with and without the cosmological constant) were studied by Zhou [21], Jing 
et.al. [25] [26] and Zhang et.al. [21]. The Dirac field around the charged string black 
hole was studied by Shu et.al. [2H]. Neutrino QNM's of a Kerr-Newmann-de-Sitter 
black hole were calculated by Chang and Shen in [29]. 

In this work we study the massless spinors around the Born-Infeld black hole. 
Born-Infeld electrodynamics was first introduced in 1930 's to obtain a classical the- 
ory of charged particles with finite self-energy [30]. In Maxwell's theory of electro- 
dynamics, the field of a point-like charge is singular at the origin and its energy is 
infinite. In contrast, in the Born-Infeld electrodynamics, the electric field of a point 
like object which is given by E r = Q/^Jr 4 + ^ is regular at the origin and its total 
energy is finite. Born-Infeld theory has received renewed attention since it turns out 
to play an important role in string theory. It arises naturally in open superstrings 
and in D-branes [31] . The low energy effective action for an open superstring in loop 
calculations lead to Born-Infeld type actions [32]- It has also been observed that the 
Born-Infeld action arises as an effective action governing the dynamics of vector-fields 
on D-branes [33]. For a review of aspects of Born-Infeld theory in string theory see 
Gibbons [M] and Tseytlin [35J. 

If we construct black hole solutions to Einstein-Born-Infeld gravity, it can be 
observed that the non-singular nature of the electric field at the origin changes the 
structure of the space-time drastically. The singularity at the origin is dominated by 
the mass term Mj r rather than the charge term Q 2 / r 2 as in the black holes of Einstein- 
Maxwell gravity [36]. Here the possibility exists for M to be positive, negative or zero 
leading to space-like or time-like singularities in contrast to the Reisser-Nordstrom 
black hole. Overall, the physical properties of black holes change drastically due to 
the non-linear nature of electrodyanmics and is worthy of study. In this paper, we 
are focusing on the static charged black holes of the Born-Infeld electrodynamics. 
This is the non-linear generalization of the Reissner-Nordstrom black hole and is 
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characterized by charge Q, mass M and the non-linear parameter (3. The Born-Infeld 
black hole was obtained by Garcia et.al. [37J in 1984. Two years later, Demianski 
[38] also presented a solution known as EBIon which differs with the one in [37] by 
a constant. Trajectories of test particles in the static charged Born-Infeld black hole 
were discussed by Breton [39J. This black hole in isolated horizon framework was 
discussed by the same author in [ID]. Gibbons and Herdeiro [H] derived a Melvin 
Universe type solution describing a magnetic field permeating the whole Universe 
in Born-Infeld electrodynamics coupled to gravity. By the use of electric-magnetic 
duality, they also obtained Melvin electric and dyonic Universes. 

The paper is presented as follows: In section 2, the Born-Infeld black hole solutions 
are introduced. In section 3, the spinor perturbations are given. In section 4, we will 
computer the QNM's and discuss the results. Finally, the conclusion is given in 
section 5. 



2 Introduction to the Born-Infeld black hole 

In this section, we introduce the Born-Infeld black hole. The action for the Einstein 
gravity coupled to Born-Infeld electrodynamics is given by, 



S = d^X\f^g 



U +L(F) 



16nG 



Here, L(F) is a function of the field strength F^ v and for Born-Infeld electrodynamics 
may be expanded to be, 



L{F)= ^ 2 (l-^l + ^j (2) 

Here, (3 has dimensions length -2 and G length 2 . In the following sections it is assumed 
that 167rG = 1. Note that when (3 — > 00, the Lagrangian L(F) approaches the one 
for Maxwell's electrodynamics given by — F 2 . 

The static charged black hole with spherical symmetry can be obtained as, 

ds 2 = -f(r)dt 2 + f{r)~ l dr 2 + r 2 (d6 2 + Stn 2 (6)d<p 2 ) (3) 

with, 



. ^ 2M , 2(3 2 r 2 / L , Q 2 \ , 4Q 2 „ (I 1 5 Q 2 > . 



Here 2F1 is the hypergeometric function. In the limit (3 — > 00, the elliptic integral 
can be expanded to give, 

m«, = 1 - ^ + % (5) 
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resulting in the function f(r) for the Reissner-Nordstrom black hole for Maxwell's 
electrodynamics. Near the origin, the function /(r) has the behavior, 



r So 



Here, 

-5 




(7) 



Depending on the values of M, Q and (3, the function /(r) for the Born-Infeld black 
hole can have two roots, one root or none. When f{r) has two roots, the behavior of 
f(r) is similar to the Reissner-Nordstrom black hole. When it has only a single root, 
the black hole behave similar to the Schwarzschild black hole. Hence the Born-Infeld 
black hole is interesting since it possess the characteristics of the most well known 
black holes in the literature. 

Extreme black holes are possible when both f{r) = and f'(r) = 0. This condition 
yields the horizon radii as, 

rex = ^ (8) 

Hence, extreme black holes are possible only if Qj3 > 1/2. We will plot the graphs for 
/(r) in the following figure. Note that for Qf3 < 1/2 only Schwarzschild type black 
holes are possible. 



4 



f 




Figure 1. The behavior of the function f(r) when the mass is changed is plotted for 
fixed values of Q and (3. In the top graph Q — 1 and [3=1. In the bottom graph 
Q = 0.3 and = 1. 

A detailed description of the characteristics of the Born-Infeld black hole is given 
and [39]. The Hawking temperature of the black hole is given by, 



T 



1 

Ait 



— + 2(3\r + i 
r+ \ 



(Q2 + r 2 ^ 



(9) 



Here, r + is the event horizon of the black hole which is a solution of /(r) = 0. The 
zeroth and the first law of black holes in Born-Infeld electrodynamics are discussed 
in detail in [13]. Static charged black hole solution to the above action with a cosmo- 
logical constant was presented in [H] [15] and was extended to higher dimensions in 



3 Massless spin | field around spherically symmet- 
ric static black hole 

In this section, we will develop the equations for a spin 1/2 field in the background of 
a static spherically symmetric space-time. The general equation for a massless Dirac 
spinor field in curved space-time can be written as, 

7 a <(^-W = (10 ) 
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Here, are the inverse of the tetrads e£. The metric tensor of the space-time 
considered and the tetrad are related by, 

g^u = Vabeffi (11) 

Here, r) a b = (—1, 1, 1, 1) is the metric in the flat space. The 7° matrices are defined 
by, 

H-^o°) (i2) 

Here, <r a are the Pauli spinors 

••-(? o> ■*-(?;) m> 

The 7 a matrices satisfy the anti-commuting relations, 

{ 7 a ,7 b } = 2^/ (14) 
The spin connection is defined as, 

r, = -i[7 a ,7 6 ]A a6 , (15) 

where, 

Kbfi = Za e bv^ (16) 

the covariant derivative of e hv given by, 

&bu;[M d^Cbi/ r^^eba (1^) 

are the Christoffel symbols. 

We will compute the Dirac equation for a space-time which is spherically sym- 
metric and static given by, 

ds 2 = -A(r) 2 dt 2 + B(r) 2 dr 2 + R(r) 2 (d6 2 + sin 2 {9)d^ 2 ) (18) 

The tetrads for this space-time are given by, 

e « = diag{A,B,R,Rsin6) (19) 

The non-zero components of A a fe M defined in eq.(16) for the above metric are computed 

as, 

A A ' 

Acm - --g 
Am - - ¥ 



• a R ' 



B 



A 23 ^ = -cosO (20) 
Therefore, the related spin connections becomes, 

r r = o 

r R' 1 2 
re= 25 77 



r v = ^-7V + ^g-^S 3 (21) 

The Born-Infeld black hole considered in this paper has the functions A = I? -1 and 



R = r. For further simplification, we will take A = B 1 = y/(r). The function ?/> is 
redefined as, 

V> = (22) 
Then the Dirac equation simplifies to, 

*» + ^(9 + iW£f4 + ^V^£ = (23) 



■^/J 5i * \drrj r \d6 2 I rsinO dip 



Before looking for solutions for the above equation, we will make another assumption. 
The field considered here is massless and has spin |. Hence the allowed solutions to 
the Dirac equation are circularly polarized. This is discussed in jlH] and [SD] at length. 
According to [SO], the above spinors have right handed circular polarization. Hence, 
the allowable spin states satisfy the condition, 



where 

Then, eq.(24) will be satisfied by, 

$ = 



(l-i 75 )$ = (24) 

7 5 = 7 o 7 y 7 3 (25) 



<$> 2 (t,r } 9,<p) 



(26) 



Now, the Dirac equation yield two identical set of equations each coupling $ x and <& 2 . 
Hence only one set of the equations will be solved. Let us redefine the components 
$i and $2 as follows: 



e _wt (27) 
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Substituting the above in eq.(23) yields, 

/iwr n rzdF\ 1 / ' cot6 TT i dH 2 dH 2 \ 1 . , 

/W„ r dG\ 1 (cot6 TT i dHx dHA 1 . . 

[77 G + r ^)F + {— Hl ~eW + ^f)jrr (29) 

To solve the angular part of the equation the following operators are defined, 

/ cot9 i d d \ . . 

a * = -[— + 7^8-^m) (30) 

_ / cot9 i d d \ . , 

d - = -{— + 7med^ + de) (31) 

These operators were discussed in [51] at length. It was shown that the operators act 
on spin weighted spherical harmonics s Yi m as ladder operators. In particular if the 
spin s = |, then the above operators gives the following relations; 

d + (_ l Y lm ) = (l + h l Y lm (32) 

2 2^ 2 

cL(i^ m ) = -(Z + \)_iY lm (33) 

2 ^ 2 

By comparing the angular parts of the equation of the Dirac equation, one can see 
that H\ and H 2 are the spin | weighted spherical harmonics given by 



Hi =i F im (34) 
# 2 =_i ^ m (35) 

2 

Hence the angular part of the eq.(28) and eq.(29) becomes A = I + |. Note that Z 
can be written as I = (2k — l)/2 where k is a positive integer. Also —l<m<l. 
Therefore, the radial part simplifies to, 

cUG - iuG + W{r)F = 

d*F + iwF + iy (r)G = (36) 
r* is the well known "tortoise" coordinate given by, 

dn = j (37) 



The function W(r) is, 



W(r) = — (38) 
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Two new functions Z± are defined as, 

Z± = F ± G (39) 

and eq.(36) can be decoupled as, 

d\Z ± + (uj 2 - V±)Z± = (40) 
Here, V± are related to W(r) as, 

V ± = ±d*W + W 2 (41) 

The effective potential V + for the Born-Infeld black hole is plotted to show how it 
changes with charge Q and the non-linear parameter f3 in the following figures. 

0.08, ■ ■ ■ ■ n 



V(r) 




2 3 4 5 6 



Figure 2. The behavior of the effective potential V + with the charge for the Born- 
Infeld black hole. Here, M — 1, f3 — 0.2 and A = 1. The height of the potential 
decreases when the charge decreases. 
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Figure 3. The behavior of the effective potential V + with the non-linear parameter 
p. Here, M — 1, Q — 1 and A = 1. The maximum height of the potential increases as 
/3 increases. The dashed one is the potential for the Reissner-Nordstrom black hole 
with same mass and charge. 
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3.1 Comparison with other fields 

Gravitational perturbation and the scalar perturbations of the Born-Infeld black hole 
were studied by Fernando in [47] [IE]- The equations were shown to be simplified to a 
Schrodinger type equation as for the spinor 1/2 field studied here. Here we will plot 
the effective potentials for the three fields to compare the behavior. 
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Figure 4. The behavior of the effective potential V+ for fields of spin 2, 1/2 and 
0. Here, M = 1, Q = 1, (3 = 0.4. The value of / = 2 for scalar and gravitational 
potential, while A = 2 for the spinor potential. 

4 Quasi-normal modes of the Born-Infeld black 
hole 

Usually, the fundamental equation of black hole perturbations given in eq.(40) cannot 
be solved analytically. This is the case for almost all black holes in 3+1 dimensions. 
In 2+1 dimensions, there are two black hole solutions ( BTZ black hole |52j and the 
charged dilaton black hole |53j) which can be solved to give exact values of QNM's. 
In five dimensions, exact values are obtained for vector perturbations by Nunez and 
Starinets [54] . 

There are several approaches to compute QNM's in the literature. Here, a semi 
analytical technique developed by Iyer and Will [55] is followed. The method makes 
use of the WKB approximation. The formula developed by Iyer and Will was for up to 
third order. Following an approach suggested by Zaslavskii [56], Konoplya presented 
the sixth order WKB formula in [57] . We will follow the formalism presented in that 
paper. The WKB formula is given by, 

Uj2 ~Il -L 2 -L 3 -L A -L,-L & =n+ l - (42) 
f 1 ^ 2 

Here, Vq and Vq represents the maximum potential and the second derivative with 
respect to the tortoise coordinates where the maximum occurs. L 2 and L3 are given 
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in [55] and L3, L4, L5 and L§ are given in [57]. n is the overtone number. The u is 
represented as u = uj r — za;/ and presented accordingly in the tables and the graphs. 
Unless states specifically, the QNM's are computed for n — 0. 

First, the quasi normal modes for V+ are computed to see the behavior with the 
non-linear parameter f3 and graphed in the following figures. 
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Figure 5. The behavior of Re u with the non-linear parameter (3 for M — 2, Q = 1 
and A = 3. 
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Figure 6. The behavior of Im u with the non-linear parameter (3 for M — 2, Q — 1 
and A = 3 

The spinor perturbation decay faster when the non-linear parameter /3 increases. 
However, beyond a certain value, the decay becomes slower as is evident from the 
graph. 

Next, we have computed the QNM's by varying the charge of the black hole. We 
also computed the QNM's for the Reissner-Nordstrom black hole with the same mass 
and the charge. The results are given in the following table and the figures. 



Q 


^r(BI) 


wj(BI) 


0J r (RN) 


u)i (RN) 


0.1 


0.57506 


0.0963636 


0.575063 


0.0963601 


0.2 


0.57796 


0.0965668 


0.578015 


0.0965175 


0.3 


0.582785 


0.0969885 


0.583084 


0.0967732 
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0.589571 


0.0976722 


0.59052 


0.0971135 
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0.598382 


0.0986511 


0.600731 


0.0975103 
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0.609334 


0.0999679 


0.614369 


0.0979046 


0.7 


0.622627 


0.10166 


0.632506 


0.0981638 


0.8 


0.638537 


0.103767 


0.657033 


0.0979584 


0.9 


0.679767 


0.109443 


0.746085 


0.0883246 
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Figure 7. The behavior of Re uj with the charge Q for M — 1, f3 — 0.04 and A = 3. 
The dotted lines show the graph for Reissner-Nordstrom and the dark lines show the 
one for the Born-Infeld black hole. 




Figure 8. The behavior of Im uo with the charge Q for M — 1, (3 — 0.04 and A = 3. 
The dotted lines show the graph for Reissner-Nordstrom and the dark lines show the 
one for the Born-Infeld black hole. 

One can observe that the decay of the spinor field in Born-Infeld black hole is faster 
than the Reissner-Nordstrom black hole for this particular values of the parameters. 
For scalar field perturbations the decay for the Born-Infeld black hole is also faster 
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than the Riessner-Nordstrom black hole [H]. Note that, for the Reissner-Nordstrom 
black holes, Q = 1, M = 1 represents the extreme black hole. This will explain why 
the uji is approaching zero around Q = 0.9 in the table. 

The behavior of the quasi normal modes for the Born-Infeld black hole with vary- 
ing A was also studied as given in the following table. 



W_R 



A 






2 


0.198638 


0.0491501 


3 


0.299979 


0.049059 


4 


0.400892 


0.0490351 
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0.501645 


0.0490251 
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0.60232 


0.0490197 
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0.702949 


0.0490165 
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0.803551 


0.0490144 




Figure 9. The behavior of Re u with the A for M = 2, (3 = 0.04 and Q = 1 
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Figure 10. The behavior of Im u with the A for M = 2, (3 = 0.04 and Q = 1 



The real part of u shows a linear relation with A while the imaginary part decreases 
to a stable limit. 
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Next, the behavior of the Im uj with temperature of the black hole is studied. The 
plot of Im(u) vs temperature is given in the Figure 11. One can see a linear behavior 
of Im uj with the temperature. This behavior is similar to the Schwarzschild-anti-de- 
Sitter black hole studied by Horowitz and Hubeny [3]. 



. 108 
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.104 

w_I 

. 102 
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.098 

0.04 0.041 0.042 0.043 0.044 
T 



Figure 11. The behavior of ujj with the Temperature for A = 3. The mass M = 1 
and /3 = 0.04. 

The overtones of the QNM frequencies can be obtained only for A > n as stated 
in [58]. The sixth order WKB method is used for A = 8, M = 2, (3 = 0.04 and Q = 1 
for the first four overtones as follows: 



n 




UJl 


1 


0.800366 


0.147305 


2 


0.794092 


0.246376 


3 


0.784922 


0.346724 


4 


0.773144 


0.448811 



The real part of the uj decreases for higher overtones, while the imaginary part in- 
creases with the overtone number n. This feature is similar to the Shwarzschild 
black hole as well as the charged black hole in Gauss-Bonnet black hole presented by 
Konoplya in [57] . 

5 Conclusions 

We have studied the spinor perturbations of a Born-Infeld black hole. The quasi- 
normal mode frequencies are computed using the sixth order WKB method. The 
parameters of the theory, charge Q, non-linear parameter /3, A and overtone number 
n are varied to compute the uj. For all values considered, the value of the imaginary 
part of uj is negative leading to the decay of the spinor field. In comparison with the 
Reissner-Nordstrom black hole, the decay rates are faster for the Born-Infled black 
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hole. It is interesting to note that for gravitational and scalar perturbations, the fields 
of the Born-Infeld black hole also decayed faster as noted in [47] [4"8] . 

One can also conclude that the Born-Infeld black holes are stable for the domain 
of the parameters considered in this paper. However, since we have done these com- 
putations for a Born-Infeld back hole which behave like the Schwarzschild black hole 
near the origin, it is necessary to do an evaluation based on all the parameters to 
fully understand the behavior of the ujj and to make remarks on stability of the black 
hole. Perhaps one can find an analytical approach to find the stability properties of 
the black hole for all parameters considered. This may be an interesting possibility 
for future work. 

A natural extension of this work is to study the stability of Born-Infeld- AdS black 
hole. The developments in string theory on AdS/CFT duality and the relation of 
Born-Infeld theory to string theory is a motivation for this. The Reissner-Nordstrom- 
AdS black hole has been shown to be unstable for linear perturbations in [55] . 

It was noted in section 2 that the Born-Infeld black hole behaves in a similar 
fashion as the Reissner-Nordstrom for certain parameters of the theory. In particu- 
lar, the possibility of extreme black holes exists. It is well known that the extreme 
Reissner-Nordstrom black hole is supersymmetric. It can be embedded in a N=2 
supergravity theory [60] [61] . Quasinormal modes of these extreme black holes were 
studied by Onozawa et.al. [62J and was shown to have the same frequencies. It will 
be interesting to see any sign of such behavior for the Born-Infeld black hole. 
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